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The  aim  of  the  present  study  is  to  investigate  the  enhancement  of  syngas  combustion 
using  the  promising  oxygen-enrichment  technology,  with  a  particular  attention  on  optimal 
operating  conditions  in  regard  to  NOx  emissions.  For  this  purpose,  a  numerical  study  is 
conducted  on  a  syngas  counter-flow  diffusion  flame,  using  air  enriched  with  oxygen  as  the 
oxidizer.  The  oxygen  concentration  ranges  from  21%  to  30%  by  volume.  Two  syngas 
compositions  are  considered  with  H2/CO  rates  equal  to  0.25  and  4,  respectively.  Flame 
structure  is  characterized  by  solving  flamelet  equations  with  the  consideration  of  radia¬ 
tion.  The  chemical  reaction  mechanism  used  is  GRI  3.0.  Computational  results  showed  that 
oxygen  addition  increases  the  flame  temperature  and  intensifies  the  radiative  heat 
transfer.  It  also  considerably  extends  flammability  limits  allowing  stable  flames  at  high 
values  of  the  scalar  dissipation  rates  and  for  lean  syngas  composition.  NOx  formation  is 
substantially  increased  with  oxygen  increment,  and  Zeldovich  mechanism  is  found  to  be 
the  main  route  of  NOx  formation.  H2-lean  syngas  flames  produce  less  NOx  at  low  scalar 
dissipation  while  H2-rich  syngas  flames  NOx  emissions  are  low  at  high  scalar  dissipation 
rates. 

Copyright  ©  2015,  Hydrogen  Energy  Publications,  LLC.  Published  by  Elsevier  Ltd.  All  rights 

reserved. 


Introduction 

Conventional  fuels  are  expected  to  be  progressively  replaced 
by  cleaner  energy  sources,  among  which,  hydrogen-rich  fuels 


such  as  syngas,  are  due  to  become  the  most  important  alter¬ 
native  ones. 

The  use  of  syngas  is  of  increasing  interest,  particularly  for 
electricity  production  and  as  a  component  of  the  Integrated 
Gasification  Combined  Cycle  concept  (IGCC). 
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The  composition  of  syngas  is  strongly  dependent  on  the 
fuel  source  and  processing  technique  [1],  Syngas  may  have 
wide  variations  of  H2/CO  ratios.  Hence,  combustion  properties 
such  as  flammability  limits  and  flame  temperature  will  be 
accordingly  variable.  In  addition,  syngas,  as  many  alternative 
fuels,  has  lower  calorific  value  than  conventional  fuels  such  as 
methane.  According  to  McLean  et  al.  [2],  a  syngas  with  a  H2/CO 
ratio  equal  to  0.5  has  a  Low  Heating  Value  (LHV)  of  10.5  MJ/kg, 
while  methane  has  a  LHV  of  50  MJ/kg  [3]. 

Oxygen-Enriched  Combustion  (OEC)  has  been  advanced  as 
cost  effective  option  to  improve  combustion  process  (stability 
and  reactivity)  [4].  Oxygen-enriched  combustion,  also  known 
as  oxygen-enhanced  combustion,  uses  oxidizers  with  oxygen 
concentrations  higher  than  that  in  air  (21%  by  volume).  It  is 
applied  since  many  years,  in  different  industries,  such  as 
cement  and  glass  factories,  waste  incineration  and  steel 
smelting.  OEC  has  received,  by  the  past,  little  attention  from 
the  academic  combustion  community,  due  to  safety  consid¬ 
erations  and  high  costs  of  oxygen  production  [4].  Recent  in¬ 
vestigations  showed  that  increasing  oxygen  concentration  in 
the  oxidizer  induces  several  benefits  including  higher  flame 
temperatures,  rising  productivity  and  significant  reduction  in 
fuel  consumption  [5,6],  Moreover,  recent  increasingly  strict 
emissions  regulations  and  costs  reduction  of  oxygen  produc¬ 
tion  have  revived  the  interest  of  researchers  to  oxygen- 
enhanced  combustion. 

Tan  et  al.  [7]  carried  out  an  experimental  study  on  the  ef¬ 
fects  of  oxygen  enrichment  in  air  and  C02  atmospheres  in  the 
case  of  methane  flame  using  a  vertical  combustor.  The  02 
concentration  was  set  to  28%  in  both  cases.  Very  high  levels  of 
NOx  were  obtained  in  02  enriched  air  combustion  due  to 
higher  flame  temperatures  while  nearly  no  NOx  emissions 
were  observed  in  02/C02  combustion  due  to  the  absence  of  N2 
in  the  feed  air.  Cheng  et  al.  [8]  investigated  by  both  numerical 
simulation  and  optical  measurements  the  structure  of  a 
planar  oxygen-enriched  methane  counter-flow  flame.  They 
particularly  focused  their  attention  on  the  effects  of  stretch 
and  oxidizer  oxygen  concentration.  They  found  that,  for  a 
given  value  of  stretch,  flame  temperature  increased  signifi¬ 
cantly  with  oxygen  concentration.  Moreover,  they  observed 
that  oxygen  enrichment  considerably  enhanced  extinction 
limits.  Chen  and  Axelbaum  RL  [9]  conducted  an  experimental 
and  a  numerical  study  to  examine  the  effects  of  oxygen- 
enrichment  and  fuel  dilution  on  flame  extinction.  They 
considered  a  counter-flow  flame  configuration  and  four 
different  fuels:  methane,  ethane,  ethylene  and  propane.  The 
range  of  oxygen  mass  fraction  variation  was  from  0.233  to  1.0. 
Results  indicated  that  an  appropriate  adjustment  of  oxidizer 
oxygen  and  fuel  concentrations  give  improved  extinction 
characteristics.  In  other  words,  an  optimal  oxygen- 
enrichment,  combined  with  fuel  dilution,  leads  to  strong 
flames  [9],  A  numerical  study  of  CHVair  and  CHV02  counter¬ 
flow  laminar  flames  was  performed  by  Urzica  and  Gutheil  [10] 
using  a  detailed  chemical  reaction  mechanism.  The  laminar 
CHV02  flame  was  studied  under  pressures  ranging  from 
0.1  MPa  to  2  MPa.  It  was  shown  that  maximum  flame  tem¬ 
perature  increased  significantly  with  oxygen  concentration  in 
the  oxidizer.  The  same  trend  was  observed  for  radicals  CO,  O 
and  OH.  Furthermore,  comparison  between  CHVair  and  CHV 
02  flames  at  strain  rate  near  extinction  revealed  that  flame 


thickness  decreased  considerably  when  oxygen  oxidizer 
content  increased.  Global  effects  of  oxygen-enhancement  on 
radiative  heat  flux  characteristics  of  methane  non-premixed 
oxyfuel  flames  were  examined  experimentally  by  Ditaranto 
and  Oppelt  [11].  The  oxidizer  oxygen  concentrations  were 
varied  from  35%  to  70%,  and  C02  was  used  as  a  diluent.  The 
measurements  performed  along  the  flame  axis  showed  that 
the  distribution  of  radiative  heat  flux  and  the  flame  length 
became  shorter  when  the  oxygen  concentration  increased. 
Moreover,  the  peak  of  radiated  heat  flux  increased  with  oxy¬ 
gen  concentration  [11].  Furthermore,  the  intensity  of  visible 
spectra  of  these  flames  was  more  important  with  oxygen 
enrichment.  Joo  et  al.  [12]  performed  a  numerical  and  an 
experimental  study  of  the  structures  of  laminar  diffusion 
methane/oxygen  and  methane/air  flames  under  a  wide  range 
of  operating  pressures  (from  1  atm  to  60  atm).  The  authors 
observed  that,  as  the  pressure  increased,  the  visible  flame 
height  of  the  methane/oxygen  flames  decreased,  while  the 
visible  flame  height  of  methane/air  flames  remained  constant. 
The  comparison  of  soot  production,  for  different  pressures, 
showed  that  maximum  soot  fraction  produced  by  methane/ 
air  flames  was  considerably  higher  than  in  methane/oxygen 
flames  [12]. 

Yepes  et  al.  [13]  analyzed  experimentally  and  numerically 
the  laminar  burning  velocities  of  a  syngas  mixture  containing 
20%  H2,  20%  CO  and  60%  N2  using  oxygen  enriched  air  with 
enrichment  levels  varying  from  21%  to  35%  by  volume.  They 
found  that  oxygen  rise  in  the  combustion  air  increased  the 
laminar  burning  velocity  of  the  syngas  mixture,  leading  to  an 
intensification  of  the  reaction  rate. 

The  literature  on  oxygen-enhanced  diffusion  flames  is 
abundant,  but  only  few  studies  are  devoted  to  syngas  oxygen- 
enhanced  flames  characteristics. 

The  aim  of  this  study  is  to  examine  numerically  the  impact 
of  oxygen  enrichment  on  the  structure  and  emissions  of 
syngas  flame  in  counter-flow  configuration.  The  calculations 
are  conducted  over  a  large  range  of  scalar  dissipation  rates 
(from  near  equilibrium  to  near  extinction)  for  oxygen  con¬ 
centration  ranging  from  21%  to  30%  by  volume,  and  for  two 
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Fig.  1  -  Flame  configuration. 
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syngas  compositions  (H2/CO  ratio  =  0.25  and  4.0).  Special  focus 
is  paid  to  NOx  emissions. 

The  article  consists  of  four  sections.  Following  a  brief 
explanation  of  the  modeling  approach,  a  general  description 
of  the  simulation  details  is  introduced.  The  numerical  simu¬ 
lation  results  are  then  presented  with  a  discussion  in  the  last 
section.  The  conclusion  summarizes  the  findings  of  the  pre¬ 
sent  work. 


Calculation  methodology 

The  flame  configuration  considered  in  the  present  work  is 
described  in  Fig.  1.  It  is  a  counter-flow,  laminar  diffusion 
flame,  stabilized  near  the  stagnation  plane  of  two  opposing 
fuel  and  oxidant  jet  flows. 

A  one  dimensional  description  of  the  local  flame  structure 
in  mixture  fraction  space  is  given  by  flamelet  equations  with 
the  scalar  dissipation  rate  as  a  free  parameter.  Flamelet 
approach  assumes  a  one-dimensional  behavior  of  the  com¬ 
bustion  phenomena  in  the  normal  direction  to  the  flame  front. 
Flamelet  equations  are  derived  from  energy  and  species 
transport  equations  by  applying  a  coordinate  transformation 
[14].  In  these  equations,  temperature  and  species  mass  frac¬ 
tions  are  function  of  a  conserved  scalar  known  as  mixture 
fraction  Z.  In  addition,  a  parameter  called  scalar  dissipation 
rate  x  totally  represents  the  influence  of  the  flow  field  on  the 
flamelet  structure: 
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The  notation  in  Equations  (1)  and  (2)  is  as  follows:  Y„,  T,  p 
are  the  nth  species  mass  fraction,  mixture  temperature  and 
density,  respectively.  Le„  is  the  Lewis  number  of  the  nth  spe¬ 
cies  defined  as  Le  =  V(pDnmCp)  where  Dnm  is  the  multi- 
component  ordinary  diffusion  coefficient,  is  the  nth  spe¬ 
cies  reaction  rate  and  %  is  the  instantaneous  scalar  dissipation 
rate  defined  by:  %  =  2Dz(VZ-VZ).  Its  modelling  is  based  on  the 
relation  below  which  is  taken  from  the  counter-flow  geometry 
[15]: 


_  j>_  g(z) 

*  “  XsVst  g(Zst) 


(3) 


Xst  is  the  scalar  dissipation  rate  at  stoichiometry  and  <p  is  a 
factor  introduced  in  order  to  include  the  effect  of  density 
variation  [15]: 


i  ^(VpoJp  + 1) 

^  4  2^/pJp  +  l 


The  subscript  oo  means  the  oxidizer  stream.  The  function 
g(Z)  is  given  as  follows  [15]: 

g(Z)  =  exp[-2(er/c-1(2Z))2]  (5) 

where  erfc is  the  inverse  of  the  complementary  error 
function. 

The  last  term  in  Equation  (2)  is  an  optically  thin  model  for 
radiative  energy  loss  from  the  flamelet.  Here,  a  is  the  Ste- 
fan-Boltzmann  constant,  P  is  the  pressure,  X„  is  the  nth 
species  mole  fraction,  a;  are  polynomial  coefficients  for  the 
Planck  mean  absorption  coefficients  [16]  and  Tb  is  the  far-field 
(background)  temperature.  Radiative  gaseous  species  consid¬ 
ered  are  CO,  C02  and  H20. 

The  validity  of  the  optically  model  for  hydrogen  and 
hydrogen  enriched  flames  (H2/He,  CO/H2/N2,  CH4/H2/N2)  was, 
previously,  confirmed  by  Barlow  et  al.  [16]  and  Frank  et  al.  [17]. 
In  addition,  many  authors  calculated  syngas  flames  using 
optically  thin  model  (Som  et  al.  [18",  Park  et  al.  [19],  Kim  et  al. 
[20]). 

The  two  input  parameters  in  Equations  (1)  and  (2)  are 
mixture  fraction  Z  and  stoichiometric  scalar  dissipation  rate 

Xst- 


Simulation  details 

Flamelet  Equations  (1)  and  (2)  are  solved  from  near  equilib¬ 
rium  to  near  extinction  until  steady  state  is  achieved 
assuming  Le„  =  1  for  all  the  species  involved  in  the  chemical 
mechanism  using  the  PrePDF  4.0  code  [21].  Computations  are 
performed  for  two  syngas  compositions  (lean  and  rich  syngas) 
using  five  levels  of  02  in  the  oxidizer.  A  total  of  ten  compu¬ 
tations  is  carried  out  (Table  1).  The  operating  pressure  is 
maintained  at  1  atm. 

The  assumption  regarding  the  Lewis  number  is  reasonable 
in  this  flame  conditions.  Differential  diffusion  was  found  to 
play  a  significant  role  in  diluted  and  non-diluted  syngas 
flames  with  high  H2  content  in  the  fuel  (>80%)  at  small  strain 
rates  (near  equilibrium)  (Park  et  al.  [19];  Drake  and  Blint  [22]; 
Park  et  al.  [23]). 

Syngas  oxidation  chemistry  is  modeled  using  GRI  3.0 
mechanism  [24]  that  involves  53  species  and  325  reactions. 
The  oxidation  chemistries  of  CO  and  H2  which  are  the  domi¬ 
nant  components  of  most  syngas  fuels  are  well  represented  in 


(4) 
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this  mechanism.  GRI  3.0  mechanism  was  used  in  many  pre¬ 
vious  studies  of  syngas  counter-flow  flames  (Park  et  al.  [25], 
Shih  et  al.,  [26]).  Its  accuracy  was  also  assessed  in  the  case  of 
syngas  in  both  air  combustion  and  oxygen-enhanced  com¬ 
bustion  employing  experimental  data  and  predictions  of 
existing  syngas  kinetics  namely  Davis  et  al.  [27]  mechanism, 
San  Diego  (Petrova  and  Williams  [28])  mechanism  and 
(Mueller  et  al.  [29])  mechanism  on  laminar  burning  velocity 
where  reasonable  agreement  was  achieved  (Yepes  and  Amell 
[13];  Som  et  al.,  [18];  Singh  et  al.,  [30]). 

In  addition,  the  study  of  Sahu  and  Ravikrishna  [31]  in 
which  five  chemical  kinetic  mechanisms  namely,  H2/CO 
mechanism  by  Davis  et  al.  [27],  Cl  mechanism  from  Li  et  al. 
[32],  full  chemistry  from  Li  et  al.  [32],  H2/CO  mechanism  by 
Ranzi  et  al.  [33]  and  GRI  3.0  [24]  were  compared  to  predict  NOx 
routes  in  syngas  flames  highlighted  the  importance  of  GRI  3.0 
mechanism  for  prompt  NO  calculation. 


Results  and  discussion 

Model  validation 

The  accuracy  of  flamelet  model  was  demonstrated  in  previous 
study  [14,15]  where  it  was  shown  that  flame  structure  is  suf¬ 
ficiently  well  described.  Further  validation  of  this  model  is 
conducted  here  in  the  case  of  oxygen-enhancement  com¬ 
bustion  employing  the  experimental  data  of  Cheng  et  al.  [34] 
which  corresponds  to  methane  as  no  measurements  of  syn¬ 
gas  combustion  with  enriched  air  in  counter-flow  diffusion 
flame  configuration  are  currently  available.  The  experimental 
device  in  Ref.  [34]  is  an  opposed-jet  burner  formed  by  two 
axisymmetric  nozzles  of  25  mm  diameter,  separated  by  a 
distance  of  12.6  mm.  The  design  of  this  burner  was  based  on 
the  study  of  Seshadri  et  al.  [35].  The  fuel  was  a  mixture  of 
nitrogen  and  methane  and  the  oxidizer  was  oxygen  at 
different  concentrations  ranging  from  23%  to  100%  by  volume. 
The  ambient  pressure  was  set  to  1  atm.  Measurements  of  the 
species  concentration  were  made  along  the  centerline  of  the 


Fig.  2  -  Comparison  between  predictions  and 
experimental  results  of  Cheng  et  al.  [34]:  peak  flame 
variations  with  oxygen  variations. 


burner  using  the  nonintrusive  Raman  diagnostic  system.  The 
local  flame  temperature  was  determined  from  the  species 
concentration  using  the  ideal  gas  law.  Fig.  2  shows  a  com¬ 
parison  between  predictions  and  experiments  for  peak  tem¬ 
perature  variations  with  oxygen  concentration  at  scalar 
dissipation  rate  value  of  60  s_1  and  methane  concentration  of 
20%  by  volume.  Very  good  agreement  is  obtained  at  this  flame 
conditions. 

Fig.  3  illustrates  a  comparison  between  predictions  and 
experiments  for  peak  temperature  variations  with  methane 
concentration.  In  this  case,  oxidizer  was  pure  oxygen  and  the 
scalar  dissipation  rate  was  60  s-1.  Here  again,  both  experi¬ 
mental  and  numerical  results  match  well.  Moreover,  the  re¬ 
sults  show  similar  extinction  limit  which  is  reached  at  8% 
methane  by  volume. 

Effect  of  scalar  dissipation  rate  and  fuel  composition  on 
maximum  flame  temperatures 

Fig.  4  shows  the  maximum  flame  temperatures,  as  a  function 
of  scalar  dissipation  rate  with  and  without  the  consideration 
of  flame  radiation,  for  both  H2-lean  (Fig.  4a)  and  H2-rich 
(Fig.  4b)  syngas  composition,  respectively.  Oxidizer  is  air,  inlet 
temperature  is  300  K  and  ambient  pressure  is  1  atm.  It  can  be 
seen  that  without  radiation,  the  maximum  adiabatic  flame 
temperature  has  a  monotonic  trend,  for  both  compositions.  It 
decreases  with  increasing  scalar  dissipation  rate,  until 
extinction.  The  adiabatic  flame  is  blown  off  because  very  high 
scalar  dissipation  rates  cause  very  short  species  residence 
times  in  the  flame  zone  and  lead  to  incomplete  reactions. 
However,  with  thermal  radiation,  the  maximum  flame  tem¬ 
perature  exhibits  a  peak  at  an  intermediate  scalar  dissipation 
rate  (x  =  5.0  s_1  for  H2-rich  composition  and  %  ~  25.0  s_1  for 
H2-lean  composition).  This  result  is  different  to  the  case 
without  radiation,  and  this  difference  grows  to  be  very  large  as 
scalar  dissipation  rate  decreases.  At  low  scalar  dissipation 
rates,  the  flame  temperature  drops  significantly  because  ra¬ 
diation  heat  losses  increase,  which  results  in  a  lower  reaction 


Methane  concentration  (%  by  volume) 

Fig.  3  -  Comparison  between  predictions  and 
experimental  results  of  Cheng  et  al.  [34]:  peak  flame 
variations  with  methane  variations. 
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rate  and  higher  reactants  leakage  [36].  The  flame  extinction 
occurs  at  high  stretch  rate  due  to  insufficient  gas  residence 
time  for  chemical  reactions.  Radiative  heat  loss  is  mainly 
dependent  on  upon  flame  volume  and  the  thickness  of  reac¬ 
tion  zone  is  relevant  to  the  reciprocal  of  scalar  dissipation 
rate.  Consequently,  the  difference  of  maximum  flame  tem¬ 
peratures  with  and  without  radiation  becomes  small  with  the 
increase  of  scalar  dissipation  rate.  This  general  trend  has  been 
identified  in  previous  studies  [37,38].  In  addition,  radiation  is 
found  to  be  important  only  at  strain  rates  lower  than  the  in¬ 
termediate  value  [37]  and  its  effect  decreases  with  hydrogen 
enrichment  [39]. 


Effect  of  oxygen-enrichment  on  flame  structure  and  NO 
emissions 

In  this  section,  the  effect  of  oxygen-enrichment  on  tempera¬ 
ture,  flame  structure  and  NO  emissions  is  examined.  Figs.  5 
and  6  illustrate  the  profiles  of  temperature  and  species  for 
oxygen  concentrations  ranging  from  21%  to  30%  by  volume  for 
two  syngas  compositions  (H2/CO  =  4.0  and  H2/CO  =  0.25). 
Ambient  temperature  and  pressure  considered  are  300  K  and 
1  atm,  respectively.  The  scalar  dissipation  rate  is  10  s_1.  From 


(a)  H2/CO  =  0.25 


(b)  hz/co  =  4,0 


x  [s'1] 

Fig.  4  -  Maximum  flame  temperature  variations  with 
dissipation  rate,  with  and  without  the  consideration  of 
flame  radiation:  (a)  lean  syngas  composition  case  (b)  rich 
syngas  composition  case. 


Fig.  5,  oxygen  addition  in  the  oxidizer  induces  a  significant 
increase  of  flame  temperature.  Indeed,  for  hydrogen  rich 
composition,  temperature  rises  by  18%  when  oxygen  con¬ 
centration  increases  from  21%  to  30%  (by  volume)  while  the 
increase  is  about  10%  for  hydrogen  lean  case.  Oxygen- 
enrichment  implies  less  nitrogen  in  the  oxidizer;  therefore, 
the  flame  dilution  is  low  in  comparison  to  air  combustion. 
Then,  heat  losses  in  flue  gases  diminish  and  less  fuel  is  needed 
to  maintain  a  given  temperature.  This  finding  signifies  that  02 
addition  in  industrial  combustion  devices  involves  an  impor¬ 
tant  fuel  saving  [40],  Moreover,  oxygen  addition  in  the  oxidizer 
induces  a  slight  shift  of  the  maximum  flame  temperature  to¬ 
wards  the  fuel  side  together  with  a  reduction  of  the  flame 
zone. 

From  Fig.  6,  oxygen  enrichment  has  an  important  influence 
on  the  syngas  combustion  kinetics.  Oxygen  concentration 
increment  intensifies  species  production  such  as  H20  (Fig.  6a), 
C02  (Fig.  6b),  OH  (Fig.  6c)  and  O  (Fig.  6d).  Fig.  6b  indicates  that 
oxygen  rise  promotes  C02  formation.  Indeed,  for  H2/CO  =  4.0 
(i.e.  rich  syngas  composition),  C02  increases  by  34%  between 
21%  and  30%  of  oxygen  concentration,  while  it  increases  by 
27%  for  H2/CO  =  0.25  (i.e.  lean  syngas  composition).  This 
feature  can  be  profitably  exploited  in  industry  to  improve  C02 
capture  in  the  flue  gas  [41]. 

Flame  fractional  temperature  is  introduced  to  characterize 
flame  radiation  as  follows  25]: 


/  = 


f  nomad  lYad 

f nomad 


(6) 


In  this  equation,  Trad  is  the  flame  temperature  with  the 
inclusion  of  radiative  heat  flux  and  Tnonrad  is  the  flame  tem¬ 
perature  without  radiative  heat  flux. 

From  Fig.  7,  the  fractional  flame  temperature  at  a  fixed 
scalar  dissipation  rate  of  10  s_1  increases  with  oxygen 
enhancement.  This  is  mainly  related  to  carbon  dioxide  pro¬ 
duction  increase.  Carbon  dioxide  being  a  strong  gas  radiation 
source.  Oxygen  addition  intensifies  radiative  heat  transfer  of 
the  syngas  flame.  This  finding  is  concordant  with  the  experi¬ 
mental  observations  of  Wu  et  al.  [40].  Furthermore,  the  effect 


Fig.  5  -  Effects  of  oxygen  enrichment  on  temperature 
profiles  (1  atm,  300  K,  x  =  10  s  1)  for  rich  and  lean  syngas. 
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Fig.  6  -  Effects  of  oxygen  enrichment  on  flame  structure  (1  atm,  300  K,  x  =  10  s  *)  for  rich  and  lean  syngas:  profiles  of  (a)  H20, 
(b)  C02,  (c)  OH  and,  (d)  O. 


of  oxygen  increment  on  the  fractional  flame  temperature  is 
amplified  by  the  increase  of  CO  mole  fraction  in  synthetic  gas 
mixture  as  CO  and  C02  are  very  radiative  species  [25]. 

Fig.  8  depicts  the  variations  of  maximum  flame  tempera¬ 
ture  with  scalar  dissipation  rate,  for  different  oxygen 


Z 

Fig.  7  -  Effects  of  oxygen  enrichment  on  radiation  factor 
(1  atm,  300  K,  x  =  10  s  1). 


concentrations.  Flame  temperature  increases  from  a  mini¬ 
mum  value,  which  corresponds  to  the  quenching  limit,  up  to  a 
peak,  and  then  drops  with  scalar  dissipation  rate  increase 
until  blow-off.  For  low  scalar  dissipation  rates,  radiative  losses 
prevail,  so  the  temperature  is  too  low  to  sustain  the  flame, 


Fig.  8  -  Effects  of  oxygen  enrichment  on  peak  flame 
temperature  as  a  function  of  scalar  dissipation  (1  atm, 
300  K). 
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leading  to  quenching  [26].  On  other  side,  the  temperature 
reduction  at  higher  scalar  dissipation  rates  is  the  result  of 
insufficient  gas  residence  time.  It  can  be  observed  that  blow- 
off  limits  are  extended  with  oxygen  enrichment,  while 
quenching  limits  remain  the  same.  For  instance,  for  H2-  lean 
syngas,  extinction  occurs  for  a  scalar  dissipation  equal  to 
2030  s_1  for  air  and  3402  s_1  and  4963  s_1  for  25%  and  30% 
oxygen  concentrations,  respectively.  In  H2-  rich  syngas  case, 
extinction  occurs  at  %  =  1296  s_1  when  oxidizer  is  air.  The 
blow-off  limits  for  25%  and  30%  of  oxygen  concentrations  are 
2570  s_1  and  4730  s_1,  respectively.  This  behavior  clearly 
demonstrates  that  oxygen  enrichment  enhances  the  stability 
of  the  syngas  flame.  These  findings  are  in  good  agreement 
with  the  results  of  Yepes  and  Amell  AA  [13]  and  Merlo  et  al. 

[41]. 

The  maximum  flame  temperature  exhibits  a  peak  at  an 
intermediate  scalar  dissipation  rate  for  a  given  value  of  oxy¬ 
gen  for  the  syngas  compositions  considered.  For  values  of 
scalar  dissipation  rate  lower  than  the  intermediate  value, 
flame  structure  is  influenced  by  combined  effects  of  adiabatic 
temperature  and  radiation  heat  loss,  whereas  only  adiabatic 
temperature  effect  prevails  at  higher  values  of  strain  rate. 

For  both  syngas  considered,  an  increase  of  the  maximum 
flame  temperature  with  oxygen  enrichment  is  noticed  at 
scalar  dissipation  rate  values  lower  and  higher  than  the  in¬ 
termediate  value.  Oxygen  addition  induces  an  increase  in 
adiabatic  flame  temperature.  This  increase  is  more  important 
than  radiation  losses  (Fig.  8)  at  strain  rates  values  lower  than 
the  intermediate  value  [13]. 

For  each  oxygen  composition,  the  maximum  flame  tem¬ 
peratures  increase  more  the  syngas  is  H2-rich  for  strain  rates 
values  below  the  intermediate  value  as  radiation  decreases 
with  hydrogen  addition.  The  opposite  behavior  is  noticed  for 
strain  rate  values  higher  than  the  intermediate  value  due  to 
the  decrease  of  adiabatic  flame  temperature  with  hydrogen 
enrichment  [26,38,39]. 

Fig.  9  illustrates  the  predicted  NO  profiles  for  the  two  syn¬ 
gas  at  different  oxygen  concentrations.  Inlet  temperature  is 
300  K,  ambient  pressure  is  1  atm  and  the  scalar  dissipation 
rate  is  taken  10  s_1.  It  can  be  observed  that  NO  emissions  in¬ 
crease  with  oxygen  and  this  trend  is  amplified  by  hydrogen 
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addition  in  the  fuel  (rich  syngas).  Moreover,  the  maximum 
values  of  NO  line  up  with  maximum  temperatures.  This  result 
is  concordant  with  the  study  of  Shih  and  Hsu  [42]. 

An  analysis  of  NO  production  pathways  (thermal  NO, 
prompt  NO  and  N20  path)  is  provided  in  Fig.  10.  NNH  route  is 
not  considered  in  this  analysis.  Indeed,  Shih  and  Hsu  [42] 
observed  that  NNH  route  becomes  the  dominant  route  for 
syngas  flames  only  at  very  low  pressures  (near  0.01  atm).  The 
results  of  Fig.  10  indicate  that  for  H2-lean  syngas,  thermal  NO 
represents  about  66%  of  the  total  NO  production  while  it  in¬ 
creases  to  about  84%  for  H2-rich  syngas.  This  signifies  that  the 
thermal  route,  also  known  as  Zeldovich's  route,  is  the  domi¬ 
nant  mechanism  for  NO  formation.  This  result  is  in  good 
agreement  with  the  findings  of  [42].  In  addition,  oxygen 
increment  is  found  to  largely  increase  the  thermal  NO 
amount.  As  a  result,  the  total  NO  production  is  substantially 
increased,  despite  the  smaller  amount  of  initial  N2  available  in 
the  oxidizer.  This  trend  was  also  noticed  by  Sung  and  Law  [43]. 

Fig.  11  illustrates  NO  Emission  Index  (EINO)  as  a  function  of 
scalar  dissipation  rate,  for  different  levels  of  oxygen  enrich¬ 
ment.  The  emission  index  EINO  is  a  global  parameter  that  has 
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Fig.  9  -  Effects  of  oxygen  enrichment  on  NO  emissions  Fig.  10  -  NO  routes  as  a  function  of  oxygen  concentration 

(1  atm,  300  K,  x  =  10  s'1).  for  (a)  a  lean  syngas  (b)  a  rich  syngas. 
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been  commonly  used  to  characterize  NO  emission  from  a 
flame.  It  is  defined  as  [44]: 

/  MNo(aNOdz 

El  NO  =  —z - ^ -  (7) 

/  (MH2<aH2dz  +  Mco^codz) 

Jo 

where  M  represents  the  molecular  weight  and  w  is  the  net 
production/consumption  rate.  EINO  exhibits  a  non-monotonic 
variation  with  scalar  dissipation  rate.  It  increases  to  a 
maximum  and  then  decreases.  The  peak  of  EINO  rises 
dramatically  with  oxygen  addition  and  this  trend  becomes 
more  noticeable  for  hydrogen-rich  syngas.  An  important 
observation  from  these  results  is  the  existence  of  operating 
conditions  that  yields  lowest  amount  of  NO  production.  In 
practice,  industrial  combustion  systems  should  operate  at  low 
scalar  dissipations  for  H2-lean  syngas  (about  5  s_1)  and  high 
scalar  dissipations  for  H2-rich  syngas  (between  100  and 
1000  s_1)  in  order  to  minimize  NO  emissions.  These  findings 
agree  with  trends  reported  in  the  literature  [25]. 


Conclusions 

The  present  study  is  a  numerical  investigation  of  the  effect  of 
oxygen  enrichment  on  the  structure  and  NOx  emissions  of  a 
counter-flow  syngas  flame  over  a  wide  range  of  operating 
conditions.  Oxygen  concentrations  vary  in  the  range  of  21%- 
30%.  Two  representative  syngas  compositions  are  considered. 
In  addition,  this  contribution  performs  an  investigation  of 
optimum  operating  conditions  for  syngas  oxygen  enriched 
combustion  in  regards  to  NOx  emissions. 

The  most  important  findings  are  the  following: 

-  Flame  thickness  decreases  with  oxygen  addition. 

-  C02  production  increase  with  oxygen-enhancement. 

-  Flame  flammability  limits  are  enhanced  with  oxygen. 

-  The  maximum  flame  temperature  shows  a  peak  at  an  in¬ 
termediate  scalar  dissipation  rate  for  a  given  value  of  ox¬ 
ygen  for  the  two  syngas  compositions  considered.  For 


Fig.  11  -  NO  emission  index  as  a  function  of  scalar 
dissipation  rate,  for  different  levels  of  oxygen  enrichment. 


values  of  scalar  dissipation  rate  lower  than  the  interme¬ 
diate  value,  flame  structure  is  influenced  by  combined  ef¬ 
fects  of  adiabatic  temperature  and  radiation  heat  loss, 
whereas  only  adiabatic  temperature  effect  prevails  at 
higher  values  of  strain  rate. 

-  For  both  syngas  considered,  an  increase  of  the  maximum 
flame  temperature  with  oxygen  enrichment  is  noticed  at 
scalar  dissipation  rate  values  lower  and  higher  than  the 
intermediate  value.  Oxygen  addition  provokes  an  increase 
in  adiabatic  flame  temperature.  This  increase  is  more 
important  than  radiation  losses  at  strain  rates  values  lower 
than  the  intermediate  value. 

-  For  each  oxygen  composition,  the  maximum  flame  tem¬ 
peratures  increase  more  the  syngas  is  H2-rich  for  strain 
rates  values  below  the  intermediate  value  as  radiation 
decreases  with  hydrogen  addition.  The  opposite  behavior 
is  noticed  for  strain  rate  values  higher  than  the  interme¬ 
diate  value  due  to  the  decrease  of  adiabatic  flame  tem¬ 
perature  with  hydrogen  enrichment. 

-  NO  emissions  dramatically  increase  with  oxygen. 

-  NO  formation  is  dominated  by  Zeldovich  mechanism. 

-  Practical  burners  like  gas  turbines  or  boilers  should  operate 
at  low  scalar  dissipations  for  H2-lean  syngas  (about  5  s_1) 
and  high  scalar  dissipations  for  H2-rich  syngas  (between 
100  and  1000  s_1)  in  order  to  minimize  NOx  emissions. 
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